Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 

Recent BaBar Studies of Bottomonium States 

V. Ziegler 

SLAC National Accelerator Laboratory, Stanford, CA, USA 



We present studies in bottomonium spectroscopy carried out with samples of (98±1) million T(2S) and (121 ±1) 
million T(3S) events recorded with the BaBar detector at the PEP-II asymmetric-energy e+e — collider at 
SLAC. Precise measurements of the branching fractions for Xbl 2(1-P) 2P) — > 7T(15) and Xw,2(2f ) — > -fT(2S) 
transitions, and searches for radiative decay to the jjf,(lS) and rn,(2S) states are performed using photons that 
have converted into an e + e~ pair. In addition, we summarize the results of a search for the spin-singlet partner of 
the X6j(lP) triplet, the hf,(lP) state of bottomonium, in the transitions T(35) — Y TT°h), and T(3S") — > 7r+7r~ hi,. 



1. Introduction 

Bottomonia are the heaviest of the qq bound states. Decays within the bottomonium family of states occur 
via 7T°, rj, u> or di-pion emission, or by electric dipole transitions or magnetic dipole transitions. Electromagnetic 
transitions between the energy levels of the bottomonium spectrum can be calculated in the quark model and are 
an important tool in understanding the internal structure of bb bound states. In particular, the measurement of 
the hyperfine mass splittings between triplet and singlet states is paramount to understanding the role of spin- 
spin interactions in quarkonium models and in testing QCD calculations. In the non-relativistic approximation, 
the hyperfine splitting is proportional to the square of the wave function at the origin, which is expected to be 
non-zero only for L = 0, where L is the orbital angular momentum quantum number of the qq system. 

In particular, the mass splitting between the T(1S*) and the rjb(lS) is a key ingredient in many theoretical 
calculations. The value measured by BaBar, m(Y(lS)) — m(r]b(lS)) = 69.3±2.8 MeV/c 2 and subsequently 
by CLEO [|| is larger than most predictions based on potential models Q, but in reasonable agreement with 
predictions from lattice calculations Q. According to ref. @, the shift in measured rib(lS) compared to QCD 
predictions might be explained by the mixing of the rjb with a CP-odd Higgs scalar. To test this model, a 
measurement of the rjb width is essential, and this was one of the BaBar motivations for searching for the rjb 
using photon conversions as described in Sec. 2. 

For L = 1, the splitting between the spin-singlet ( 1 Pi) and the spin- averaged triplet state (( 3 Pj)) is expected 
to be AMhf = M{ 3 Pj) — M(}Pi) ~ 0. The 1 P\ state of bottomonium, the hb(lP), is the axial vector partner 
of the P-wave Xbj{lP) states. Its expected mass, computed as the spin-weighted center of gravity of the 
XbjO-P) states Q, is 9899. 87± 0.27 MeV/c 2 . Higher-order corrections might cause a small deviation from 
this value, but a hyperfine splitting larger than 1 MeV/c 2 might be indicative of a vector component in the 
confinement potential H. The hyperfine splitting for the charmonium X P\ state h c is measured by the BES 
and CLEO experiments |9Hll| to be ~0.1 MeV/c 2 . An even smaller splitting is expected for the much heavier 
bottomonium system Q. 

The BaBar [l^] experiment collected data at the narrow T(nS) resonances in order to carry out detailed 
studies in bottomonium spectroscopy 1 . 

We summarize the results of a study of radiative bottomonium transitions using 7 — > e + e~ conversions, and 
of searches for the h b (lP) state in T(35) -> TT+n-h b (lP) decay pj] and T(35) -> 7r%(lP) decay 0]. 



2. Study of radiative bottomonium transitions using 7 — > e + e~ conversions 

The rate of 7 —5- e + e~ conversion in the detector material, and of the reconstruction of the e + e~ pairs is 
much lower than that for photons reconstructed using the BaBar electromagnetic calorimeter. However, the 
substantial improvement in energy resolution for photon conversions results in better separation of photon 
energy lines. 



1 These samples correspond to 14 fb 1 and 30 fb 1 of T(2S) and T(3S I ) data, respectively. 
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Figure 1: (left) End view of the the BaBar inner detector along the beam axis as seen by converted photons. Points 
indicate the number of converted photon vertices per cross-sectional area, as measured in a subset of the "test sample" 
data. From the center outwards, features of note include the beam pipe, the SVT [e.g. hexagonal inner layers) and its 
support structure rods, the support tube, and the inner wall of the DCH. (right) The inclusive converted photon energy 
spectrum from (a) Y(3S) and (b) T(2S) decays. The shaded areas indicate different regions of interest considered in 
detail in this analysis. 



Photon conversions are reconstructed with a dedicated fitting algorithm that pairs oppositely charged particle 
tracks to form secondary vertices away from the interaction point. To remove internal conversions and Dalitz 
decays, and to improve signal purity, the conversion vertex radius is restricted to the beampipe, SVT, support 
tube, and inner wall of the DCH (see Fig. 1). 

The data are divided into four energy ranges (see Fig. 2) corresponding to bottomonium transitions, in Y(35) 
data, 

• 180 < E*(i) < 300 MeV: x&j(2P) -> 7T(25) 

• 300 < E*(<y) < 600 Mev: T(35) -> JXbj(lP) and T(35) -> yr) b (2S) 

• 600 < E*(<y) < 1100 MeV: Xb.j{2P) -> tT(15) and T(35) -> 7%(15) 
and in T(25) data, 

• 300 < £7* (7) < 800 MeV: Xbj0-P) ~> 7T(15) and T(25) 7?y 6 (lS'). 

The photon energy (E*(j)) is analysed in the center-of-mass (CM) frame of the initial T(2S, 35) system. 
The photon spectra from subsequent boosted decays (e.g. Xbj{nP) —> y¥(lS)) are consequently affected by 
Doppler broadening due to the motion of the parent state in the CM frame. This broadening of the lineshapes 
of the various photon lines is incorporated into the fits used to extract the yields for the branching fractions 
obtained. 
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Figure 2: Fit results for the inclusive converted photon spectrum after background-subtraction for (left) Y(3S), and 
(right) T(2S) data. The purple, blue, green, gray, and cyan curves represent (left) Xt>o,i,2(2P) — ¥ yT(lS) and (right) 
Xm,i,2(1-P) — >• 7T(15'), initial state radiation (ISR), and T(nS) — > r yrjb(lS)events, respectively. 



Transition Branching Fraction Measurement (%) 



Xbo(2P) 


7T(2S) 


-4.9±2.9to;g ±0.5 (< 2.9) 


Xbi(2P) 


-> 7 T(2S) 


19.5 ± l.l±l:o ± 1-9 


Xb2(2P) 


->■ 7T(2S) 


8.6±^ ±0.5 ±1.1 


T(3S) -> 


7X62 (IP) 


10.6 ±0.3 ±0.6 


T(3S) -> 


7X61 (IP) 


0.5 ±0.3l^ (< 1.1) 


T(35) -> 


7X60 (IP) 


2.7 ±0.4 ±0.2 


Xw(l-P) 


->7T(1S) 


2.3 ± 1.5±J;? ± 0.2 (< 4.6) 




^ 7 T(1S) 


36.2 ±0.8 ±1.7 ±2.1 


X62(1P) 


->7T(1S) 


20.2 ± 0.7±i;° ± 1.0 


T(25) - 


>7»?6(15) 


0.11 ± 0.04lg;^ (<0.22) 


Xw(2P) 


-^ 7 T(1S) 


0.7±0.4±o.i ±0.1 (< 1.2) 


Xw(2P) 


^ 7 T(1S) 


9.9 ±0.3 ±0.4 ±0.9 


Xm(2P) 


->7T(1S) 


7.1 ±0.2 ±0.3 ±0.9 


T(3S') - 


>7%(1S) 


0.059 ± 0.016j;g ^e 



Table I: Summary of Branching Fraction Measurements. The first uncertainty is statistical and the second systematic. 
Upper limits are given at the 90% confidence level. 



The number of signal events for a given bottomonium transition is extracted from the data by performing a 
X 2 fit to the E*(pf) distribution in 1 MeV intervals. The functional form and parameterization for each photon 
signal are determined from simulation. In general, the lineshape is related to a Crystal Ball function [26| . i.e. a 
Gaussian function with a power- law tail. This functional form is used to account for bremsstrahlung losses of the 
e + e _ pair. The underlying smooth inclusive photon background is described by the product of a polynominal 
and an exponential function. This functional form adequately describes the background in each separate energy 
range. 

The decays T(3S) — > 7Xbo,2(l-P) are observed, and precise measurements of the branching fractions for 
Xm,2(1-P, 2P) — > 7T(15) and x&i,2(2P) — > 7Y(2S') decays obtained. The background-subtracted fit results are 
shown in Fig. 2. The product of branching fractions B(T(nS) — > 7Xfcj) x B(Xb.J —> T(15)) obtained from the 
fit of Fig. 2 are consistent with, and improve upon, the current values 

The searches for rjb(lS) and r)b(2S) states using the converted photon energy spectrum are largely inconclusive. 
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Over a range of approximately 9974 < m Vb(2 s) < 10015 MeV/c 2 , we find S(T(3S") 7%(25)) < f.9 x 10~ 3 _[at 
90% C.L.). This value is consistent with, but does not improve upon, the upper limit obtained by CLEO [15]. 
Due to low efficiency and high background levels, no evidence for T(25) — > jr)b(lS) is found. The most 
significant peaking structure seen in the E*(^y) energy region expected for the T(3S) — > jrib(lS) transition, if 
interpreted as an rjb signal, trends toward recent potential model [l6[ and lattice [T3] predictions. However, the 
small (< 3<r) signicance of the result is insufficient to measure the rjb mass from the present analysis. Taking 
advantage of the improved resolution from a converted photon technique to make a definitive measurement of 
the rjb mass and width will require much more data from future experiments. 

The measured branching fractions obtained from the fits of Fig. 2 are listed in Table 1. 



3. The BaBar search for the h b (lP) in T(35) data 

The hb(lP) state is expected to be produced in T(35) decay via tt° or di-pion emission, and to undergo a 
subsequent El transition to the r) b (lS), with branching fraction (BF) B(h b (lP) -> 7%[l.S]i ~ (40 - 50)% 0. 
The isospin- violating decay T(35) — > n°hb{lP) is expected to have a BF of about 0.1% |18l [l9i . while theoretical 
predictions for the transition T(35) -> n + w~h b (lP) range from ~ 10" 4 [3 up to ~ 10~ 3 |20l-[22T]. 

We search for a signal in the inclusive recoil mass distribution against di-pion {ran) or 7r° (m leco ii(ir )) 
candidates. 

3.1. The search for the h b (lP) in the decay Y(3S) — > n + n-h b 

The hb signal is expected to appear as a peak in the m# distribution on top of a smooth non-peaking 
background from continuum events (e + e~ — > qq with q — u,d,s,c) and bottomonium decays. Several other 
processes produce peaks in the recoil mass spectrum close to the signal region. Hadronic transitions T(35) — ) 
tt + tt~T(2S) (hereafter denoted T 3 ^ 2 ) produce a peak centered at the T(25) mass (m[Y(25)] = 10.02326 ± 
0.00031 GeV/c 2 0). The cascade process T(35) -4 T(2S)X, T{2S) -> tt+tt-T(1S) (T 2 ^ 1 ) results in a peak 
centered at 9.791 GcV/c 2 . The peak is offset from the T(IS') mass by approximately the T(3S I ) to T(25) mass 
difference. Doppler shift and broadening further affect the position and width of this peak. When the T(25) 
parent in T 2 ^ 1 decays is produced through the T 3 ^ 2 channel, a pion from the T(3S I ) decay can be combined 
with an oppositely-charged track from the T(25 l ) decay to produce a broad distribution centered around 9.9 
GeV/c 2 . The T(2S) is also produced through the initial-state radiation (ISR) process e + e~ — > jisrT(2S) 

(T 2 ^jj). Of the nine possible T(35) — > XbJii^P)l, Xbj'i^P) — > tt + it~ XbjO-P) decay chains {xl J )> om y th° se 
for J' = J = {1,2} have been reported Q; these should generate two slightly separated peaks near 9.993 
GeV/c 2 , while the contributions with J' ^ J or with J = are expected to be negligible. 

No hb(lP) signal is seen in the T(35) — > n + 'K~hb production channel. Assuming the hb mass to be 9.900 
GeV/c 2 , we set a 90% C.L. upper limit By < 1.2 x 10~ 4 . We exclude, at 90% CL, values of By above 1.8 x 10~ 4 
for a wide range of assumed hb mass values. These results disfavor the calculations of Refs. f2(ll-4"2^| . Similarly, 
a recent measurement of the T(l 3 Dr) — > T(lS)n + ir~ branching fraction [23| disfavors the calculations of 
Ref. [HI [22| ■ The predictions of Ref. [24[ are at least one order of magnitude smaller and are not contradicted 
by our result. 

Figure 3 shows the distribution of ma after subtraction of the non-peaking background. An expanded view 
of the Xb J region is shown on the right. 

A branching fraction measurement B[T(35) 7r+7r"T(2S')] = (3.00 ± 0.02(stat.) ± 0.14(syst.))% improves 
over the current world average (2.45 ± 0.23)% 0- Product branching fractions S[T(35) ->■ Ax&i(2P)] x 
B[ X bi(2P) -> n+ir-xbi] = (1.16 ± 0.07 ± 0.12) x KT 3 , B[T(35) A X62 (2F)] x B[x 62 (2P) -»• vr+^xd = 
(0.64 ± 0.05 ± 0.08) x lO" 3 , and B[T(35) -> XT(2S)} x B[T(2S) -> tt+tt-T] = (1.78 ± 0.02 ± 0.11)% are also 
obtained from this analysis. In addition, we obtain a measurement of 331.50 ± 0.02(stat.) ± 0.13(syst.) MeV/c 2 
for the T(35 , )-T(25) mass difference. 

3.2. Evidence for the h b (lP) in the decay T(35) ir°h b 

Evidence for the hb(lP) state is obtained in the decay T(35) — > TT°hb(lP), by requiring a photon with an 
energy consistent with that for the hb(lP) — > jr]b(lS) transition. The number of ir° events in m roco ii(7r°) is 
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determined from a fit to the to 77 distribution in each ?7i roC oii(7r ) interval. 

A 7T° candidate is reconstructed as a photon pair with invariant mass 771(77) m the range 55-200 MeV/c 2 (see 
Fig. 4). In the calculation of rn reco n(ir°) , the 7-pair invariant mass is constrained to the nominal ir° value [7[ in 
order to improve the momentum resolution of the ir° . We employ a simple set of criteria based on kinematic and 
event shape variables to suppress backgrounds. Photons from ir° decays are a primary source of background in 
the region of the signal photon line from h\, —> jrjb transitions. Similarly, many misreconstructed tt° candidates 
result from the pairing of photons from different 7r°'s. We impose the ir° veto 2 condition on all photon candidates 
to improve the purity of the ir° and photon candidate sample. 

We obtain the m reC oii(7r ) distribution in 90 intervals of 3 MeV/c 2 from 9.73 to 10 GeV/c 2 . For each m reco ii(7r ) 
interval, the 771.(77) spectrum consists of a ir° signal above combinatorial background (see Fig. 4). We construct 
the 77i re coii(7r ) spectrum by extracting the ir° signal yield in each interval of TrirecoiK 71 " ) from a fit to the 
777,(77) distribution in that interval. The 77i r ecoii("" ) distribution is thus obtained as the fitted ir° yield and its 
uncertainty for each interval of m roco ii(7r ). 

We use the MC background and MC 7r°-signal distributions directly in fitting the 777(77) distributions in data. 
For each 777 re coii("" ) interval in MC, we obtain histograms in 0.1 MeV/c 2 intervals of 777(77) corresponding to 
the 7r°-signal and background distributions. The 7r°-signal distribution is obtained by requiring matching of 
the reconstructed to the generated 7r°'s on a candidate-by-candidate basis (termed "truth-matching" in the 
following discussion) . The histogram representing background is obtained by subtraction of the 7r° signal from 
the total distribution. 

For both signal and background the qualitative changes in shape over the full range of m leco ii(Tr°) are quite 
well reproduced by the MC. However, the tt° signal distribution in data is slightly broader than in MC, and is 
peaked at a slightly higher mass value. The 777(77) background shape also differs between data and MC. To 
address these differences, the MC ir° signal is displaced in mass and smeared by a double Gaussian function 
with different mean and width values; the MC background distribution is weighted according to a polynomial 
in 777(77). The signal-shape and background-weighting parameter values are obtained from a fit to the 777(77) 
distribution in data for the full range of 777 reco ii(7r°). At each step in the fitting procedure, the ir° signal and 
background distributions are normalized to unit area, and a \ 2 between a linear combination of these MC 
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Figure 3: (left) Results of the search for T(35) —5- tv + tt~ hb(lP). The ran spectrum after subtraction of the continuum 
background component. The curves represent the T\sr (dotted), K% (double-dot-dashed), xt' 1 (dashed), and Xb' 2 (dot- 
dashed) components. Inset: expanded view in the hb region after subtraction of continuum and peaking backgrounds. 

(right) The ma spectrum in the xt J re gi° n after subtraction of continuum and Ks background components: points 



represent data, while the curves represent the fitted model (solid), the xt' 1 (dashed), xt' 2 (dot-dashed), and T J 
(long-dashed) components. 



2 A signal photon candidate is rejected if, when combined with another photon in the event (72), the resulting 772 invariant 
mass is within 15 MeV/c 2 of the nominal tt° mass; this is called a it veto. 



G 
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Figure 4: The result of the fit to the 771(77) distribution in data (data points) for the full range of m reco ii(7r°). The solid 
histogram shows the fit result, and is essentially indistinguishable from the data; the shaded histogram corresponds to 
the background distribution. 



histograms and the 771(77) distribution in data is computed. The fit function provides an excellent description 
of the data (x 2 /iVDi v =1446/1433; A r _DF=Number of Degrees of Freedom) and the fit result is essentially 
indistinguishable from the data histogram. The background distribution exhibits a small peak at the ir° mass, 
due to interactions in the detector material of the type nir + — > pir° or pn~ —¥ nir° that cannot be truth-matched. 
The normalization of this background to the non-peaking background is obtained from the MC simulation, which 
incorporates the results of detailed studies of interactions in the detector material performed using data [2a |. 
This peak is displaced and smeared as for the primary 7r° signal. 

The fits to the individual 771(77) distributions are performed with the smearing and weighting parameters 
fixed to the values obtained from the fit shown in Fig. 4. In this process, the MC signal and background 
distributions for each 77i reC oii(7r ) interval arc shifted, smeared, and weighted using the fixed parameter values, 
and then normalized to unit area. Thus, only the signal and background yields are free parameters in each 
fit. The \ 2 n t to the data then gives the value and the uncertainty of the number of ir° events in each m reco ii 
interval. 

To search for an hb signal, we perform a binned x 2 fit to the m reco n(ir°) distribution obtained in data. The 
/if, signal function is represented by the sum of two Crystal Ball functions [26[ with parameter values, other 
than the hb mass, m, and the normalization, determined from simulated signal T(3S) —t ir°hb events. The 
background function is a fifth order polynomial, where parameter values are obtained from a fit which excludes 
the signal region (see Fig. 5(a)). 

We then perform a fit over the signal region to search for an hb signal of the expected shape, taking account 
of the correlated uncertainties related to the polynomial interpolation procedure, but with the fitted values of 
the polynomial coefficient fixed. 

In Fig. 5(b) we plot the difference between the distribution of m reco ii(7r°) and the fitted histogram of Fig. 5(a) 
over the entire region from 9.73 GeV/c 2 to 10.00 GeV/c 2 ; we have combined pairs of 3 MeV/c 2 intervals from 
Fig. 5(a) for clarity. 

The width of the observed signal is consistent with experimental resolution, and its significance is 3.3 a, 
including systematic uncertainties. The measured mass value, m = 9902 ± 4(stat.)±2(syst.) MeV/c 2 , is 
consistent with the expectation for the hb(lP) bottomonium state, the axial vector partner of the x&j(1-P) 
triplet of states. We obtain the value (4.3 ± l.l(stat.) ± 0.9(syst.)) x 10 -4 for the product branching fraction 
B(T(3S) -> n h b ) x B(h b -> 1Vb ). 

This measurement combined with the expected rate for the hb(lP) — > jrib^S) transition Q would indicate 
that the T(35) — > Tr + TT~hb(lP) is suppressed by a factor greater than 3 w.r.t. T(3S) — > n°hb(lP). 
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Figure 5: (a) The m rcco ii(7r°) distribution in the region 9.73 < m rcco ii(7r°) < 10 GeV/c 2 for data (points); the solid 
histogram represents the fit function described in the text. The data in the ht signal region have been excluded from the 
fit and the plot, (b) The m rcco ii(7r°) spectrum after subtracting background; in the hb signal region the data points are 
shown as squares, and the area with diagonal shading represents the uncertainties from the background fit; the shaded 
histogram represents the signal function resulting from the fit to the data in the signal region. 



4. Conclusions 



A study of radiative transitions between bottomonium states using photons that have been converted to e + e~ 
pairs in the detector material led to the observation of T(35) — > 7Xfe0,2(l-P) decays, and to precise measurements 
of the branching fractions for Xbi,2^-P, 2P) — > jT(1S) and Xfci,2(2P) — > jT(2S) decays. 

Studies of inclusive di-pion and ir° transitions was carried out using T(3S f ) events. While we find no evidence 
for the bottomonium spin-singlet state /i&(1P) in the invariant mass distribution recoiling against the tt + it~ 
system, we obtain evidence at the 3.3c level for the hb(lP) in T(3S f ) — > Tr°hb(lP) decay. Shortly after the BaBar 
hi, search results were obtained, preliminary results of a search for the hb in the reaction e + e~ —> hb(nP)-n + n~ 
(n = 1,2) in data collected near the T(5S) resonance were announced by the Belle Collaboration [27|. The 
hb(\P) mass measured therein agrees very well with the value reported in this paper. It is interesting to note 
that the e + e~ — > hb(nP)ir + ir~ cross sections measured in the Belle analysis relative to the e + e~ — > T(25')7r + 7r _ 
cross section are large, indicating an anomalous production rate for the hb(lP) and hb(2P), and suggesting that 
these states are maybe produced via an exotic process that violates the suppression of heavy quark spin-flip. 
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